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SOFTWARE DEVELOPMENT FOR INFRASOUND 
MEASUREMENT SYSTEM 


By 

Camille Khalaf 1 and John W. Stoughton 2 


I. INTRODUCTION 

A. General Description 

This manual is a brief description of mainly a software package devel- 
oped for detection and analysis of infrasounds produced by air turbulence. 
The software package operates on the APPLE®* computer with its peripherals, 
while infrasounds (.1 Hz to 10 Hz) are captured by a condenser microphone 
system deploying an array of 7 microphones developed by the NASA Acoustic 
and Vibration Instrumentation Section (AVIS). 

The aim of this project is to identify and locate severe weather storms 
by infrasound as produced by air turbulence. Successful detection and 
location of severe weather storms will provide a positive basis for the de- 
tection and location of large scale clear air turbulence. The signal 
analysis techniques to be employed are based on the results provided by our 
software package DAISE (Digital Analysis of Infrasound Experiments) . For 
more information about the microphone system and the auto-correlation 
technique, please refer to the NASA publication NASA TN D-8327 [1]. 

The following is a summary of what DAISE provides: 

1. Samples raw data from the microphone system with a variable 
sampling frequency and stores it permanently on floppy disks. 

2. Displays raw data, from any single microphone, on the video screen 
using the high resolution graphics (HGR) mode. 


1 Graduate Research Assistant, Department of Electrical Engineering, Old 
Dominion University, Norfolk, Virginia 23508. 

2 Associate Professor, Department of Electr’.cal Engineering, Old Dominion 
University, Norfolk, Virginia 23508. 

*This trade name is used for descriptive purposes only and the authors do 
not intend any endorsement of the product. 



3. Analyzes the data through computing: 

a. Amplitude spectrin of any channel as phase and magnitude. 

b. Power or cross-power spectrum between any 2 channels as phase 
and magnitude. 

c. Auto-correlation function of any channel or cross-correlation 
between any two. 

4. Displays any of the above functions using HGR mode . 

5. Saves results by printing any of the above displays. 

B. System Description 

a. System 3 lock Diagram 

A block diagram of the hardware system is illustrated in Figure 1. All 
parts are essential for proper operation. However, the user can do without 
the SILENTYPE®* printer if he does not wish to print any displays generated 
by the software. A description of the individual parts follows below. 

b. Microphone System 

This system provides the APPLE®, through an A/D converter, with the raw 
data needed for analysis. The number of outputs (microphones) could be as 
many as 8, the maximum channels the A/D can handle at one time. 

Note that as long as the rest of the APPLE® system is concerned, this 
block could be any external source providing an output ranging from 0 to 5 
volts . 

c. A/D Converter 

This block is the way of communication between the APPLE® and the out- 
side world. It is an 8 bit Analog to Digital converter, that takes as many 
as 8 analog inputs, ranging from 0 to 5 volts, and converts them to an 8 bit 
digital output (0 to 25 volts) . When specified by the user, the APPLE® 
scans the A/D and, the digital data is transferred into its internal 
memory . 

Note , for proper operation, any channel that is not used on the AID 
when sampling has to be grounded , otherwise the data on the used channels 
will be distorted. The AID unit has to be inserted in slot #4 in the back 
of the APPLE® monitor. 

*This trade name is used for descriptive purposes only and the authors do 
not intend any endorsement of the product. 
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Figure 1. System Block Diagram 









d. APPLE® Monitor 


Obviously, this is the processing unit of the whole system and it com- 
municates with all other blocks. It could be one of the following: 

APPLE® computer - 64 K 

APPLE II® computer - 64 K 

APPLE II PLUS® computer - 64 K 

If you remove the cover of your APPLE® , you will see 6 peripheral slots in 
the back, that is where all the other system blocks have to be connected. 
Please see the appropriate block description for the slot number in which it 
should b a connected. 

Note : Make sure to turn the power off before you try to insert or 

remove any of the peripheral cords. 

e. Disk Drives 

Two disk drives are required in order to run the software package. 

These 2 drives are connected through one peripheral card in one slot, 
typically slot #6. They use floppy diskettes on which programs and data are 
stored and retrieved by the APPLE®. 

i . Video Display 

This block is normally VIDEO 300 provided by APPLE® Computer Inc. , but 
it could be any comparable video system. It must be connected through slot 
# 0. The user interaction with the system takes place through the video 
screen as well as plotting results using the high resolution graphics capa- 
bility of the APPLE®. 

g. SILENTYPE® Printer 

The SILENTYPE® is the only type of printer that can be interfaced to 
the system and interact with our software, and must be connected through 
slot #1. The printer is necessary in the sense that it is the only way of 
saving your results as printed displays for later use. 

II. SOFTWARE DESCRIPTION 

A. DAISE Block Diagram 
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The software structure is very straightforwaro and can be easily 
understood by the aid of the diagram on Figure 2. As shown, DAISE consists 
of 3 levels of subroutines. The first level includes the preliminary w>rk 
th. i. has to be done before any analysis, the second level is where all the 
computations take place, while displaying the results of the computations 
is directed through level 3. Any level is reached only from the previous 
one and whenever the command EX is entered the control is transferred to the 
previous level again. Upon execution of any subroutine DAISE will return to 
the same level from which the subroutine was called. In the next 3 sections 
the subroutines of all 3 levels are fully described. 

B. Start-Up Procedure 

a. Diskettes needed 

In order to start up the system in addition to the hardware described 
in section I.B., at least three version 3.3 diskettes are required: 

1. DAISE diskette containing the following binary files: 

- PASS0 

- PASSl in addition to one text file: 

- PASS2 -MASTER 

- RUNTIME 

- ADC .0 

- DAISE. OBJ 

2. GPLE diskette, containing two binary files: 

- PLE DOS MOVER 

- PLE. DM 

3. DATA diskette that may or may not contain any files, but it has to 

be at least initialized. 

b. Step by Step Procedure 

Step One 

Insert DAISE diskette in disk drive 1 and GPLE diskette in disk drive 
2, and then boot the system. 


5 



COMMAND 


ORIGINAL PA... li 
OF POOR QUALITY 



Figure 2. Block Diagram Description of DAISE 










Step Two 

Type the command: EXEC MASTER, D1 (cr). This command will execute an 

executive file that moves DOS to the additional 16K memory card and then 
loads the compiler run time table for the APPLE* monitor. 

Step Three 

When indicator lights on both disk drives go off, remove GPLE diskette 
from disk drive 2 and insert DATA diskette instead. 

Step Four 

Type the command: 

BRUN DAISE.OBJ, D1 (cr) 

The diskette in drive 1 will spin for a few seconds, and then the first 
three lines of DAISE will be printed on the video screen leaving you with 
the first level subroutines. 

C. First Level Subroutines 
a. SAMPLE command: SA 

Description 

This subroutine samples raw data from the microphone system through an 
A/D converter. The data is an 8 K byte block sampled at 8 channels (micro- 
phones), 1024 bytes per channel. The subroutine scans all 8 channels 
simultaneously and then stores the data on the disk (Drive 2) as 8 binary 
subfiles . 

User Intervention 

The user has to intervene twice. First, right after he gives the 
command SA, he will be asked to input the sampling frequency constant. At 
this point the APPLE* will be waiting for a real nunber as input. (See log 
graph in Figure 3 for range and frequency correspondence) . Second, right 
before storing data on the disks the user is asked to input or file name as 
a string variable. Thif. variable will be common to all 8 subfile names as: 
"File"-c"x" where "File" stands for file name specified by the user, x is an 
integer in the range 1-8 and c designates a channel. 



SAMPLE FREQUENCY (Hz) 




After you load and run 9AISE, the following It dit played on your video 
screen: 

1 TYPE COMMAND SA, PR, AM, CA, EX, FOR: 

2 SAMPLE, PREVIEW, ANALYSIS, CATALOG, EXIT 

3 ? 

AT THIS POINT TYPF "SA" FROM THE KEYBOARD, THEN YOU'LL SEE: 

4 INPUT SAMPLING FREQUENCY CONSTANT 

5 ? 

NOW TYPE "68". THE APPLE* WILL SPEND 20 SECONDS SAMPLING (51.2 He) 

BEFORE REQUESTING: 

6 INPUT FILE NAME 

7 ? 

TYPE "TEST". NOW THE APPLE WILL STORE 8 SUBFILES UNDER THE FOLLOWING 
NAMES: 

TEST-C1 

TEST-C2 

TEST-C3 

TEST-C4 

TEST-C5 

TEST-C6 

TEST-C7 

TEST-C8 

Lin«s 1, 2, end 3 are printed on the video screen again, then waits for 
a new command. 

Notes 

- If an invalid cocmand is typed in line 3 above, DAISE will simply refuse 
it, flint the first 3 lines, and wait for a valid command. 

b. PREVIEW, command: PR 

Description 

The PREVIEW subroutine allowi a preview of the raw data sampled and 
stored on the disk at any earlier time. It displays the data on the video 
screen using the high resolution graphics display of the APPLE*. The 


displayed data consists of only 256 samples out of 1024 belonging to one 
subfile (channel). To preview another quarter, the use’' has to call PREVIEW 
another time. The vertical axis on the display is the voltage axis while 
the horizontal one is the time axis. 

User intervention 

At the veiy beginning of the subroutine the user is asked to input a 
subfile name. This name, as explained earlier, is of the fora "Fi le-name"-c"x M 
in which the file and the channel are specified. Then the user is asked to 
input the desired quarter as an integer in the range 1 to 4, representing the 
first, second, third and fourth 256 sample blocks of the subfile. After these 
two inputs the axis will be drawn on the display and the APPLE® will wait for a 
scaling factor for the data. The scaling factor could be any real number by 
which every data value will be divided before it is printed as a dot on the 
display. An optimum value for the scale is 1 since this is the value that 
scales the vertical axis between 0 and 5 volts. 

Example 

1 TYPE COMMAND SA, PR, AN, CA, EX, FOR: 

2 SAMPLE. PREVIEW. ANALYSIS. CATALOG.EXIT 

3 ? 

PR (cr) 

4 INPUT FILE AND CHANNEL 
? 

TEST-C1 (cr) 

5 DESIRED QUARTER? 

3 (cr) 

(THE AXIS ARE DRAWN) 

6 SCALING FACTOR? 

7 

1 (cr) 

Now the data will be displayed and th? APPLE® goes back to lines 1, 2, and 3 
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above. 


c. ANALYSIS, command: AN 

Description 

This subroutine includes most of the program since it is the route to 
all DAISE analysis features. It provides complete spectral as well as cor- 
relation analysis of the raw data. By selecting this subroutine you could 
compute and display the following: 

- Ampl. spectrum of any channel as phase and MG. 

- Power spectrum of any channel as phase and MG. 

- Cross-power between any 2 channels as phase and MG. 

- Printing real and imaginary values of any spectral comp, of ampl., 
power or x- power spectrum 

- Auto-correlation function of any channel 

- X-correlation between any 2 channels 

In addition, you could print any of the above displays on the SILENTYPE® 
printer. 

ANALYSIS leads you to the second level of the block diagram in which you 
have the choice of executing four different subroutines: SET-UP, COMPUTE, 

DISPLAY, EXIT. These subroutines will be studied later in this manual. 

User intervention 

There are many interventions that will not be listed here since they 
will be listed in second level subroutines. The major decision in ANALYSIS 
is choosing one out of the 4 subroutines listed above. You could chose any 
route at any time but for reasons that will be clear later the normal order 
of execution is: 

SET-UP, COMPUTE, DISPLAY and then EXIT, please read about these sub- 
routines before you try to execute any of them. 

Example 

1 TYPE COMMAND SA, PR, AN, CA, ER, FOR: 

2 SAMPLE, PREVIEW, ANALYSIS, CATALOG, EXIT 
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? 

AN (cr) 

4 TYPE COMMAND SE, CO, DI, EX, FOR: 

5 SET-UP, COMPUTE, DISPLAY, EXIT 

6 ? 

SE (cr) 

Notes 

An invalid command in line 6 above will be ignored, lines 4, 5 and 6 
will be printed again, and DAIS will be waiting for a valid command. 

d. CATALOG, command: CA 

Description 

This command lets you look at the contents of the disk in Drive 2, where 
your data files are stored. The object of this subroutine is to review the 
file names so you can choose the right file to preview or analyze later. 
This is done by a normal catalog command of APPLE 9 software. 

User intervention 

If the file names do not fit in one screen, just press RETURN to look 
at the rest of the files. This can be repeated until DAISE returns the 
beginning of the program. 

Example 

1 TYPE COMMAND SA, PR, AN, CA, EX, FOR: 

2 SAMPLE, PREVIEW, ANALYSIS, CATALOG, EXIT 

3 ? 

CA (cr) 

4 DISK VOLUME 254 

5 TEST-C1 

6 TEST-C2 

• • • • 

• • • • 

12 TEST-C8 
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13 TYPE COMMAND SA, PR, AN, CA, EX, FOR: 

14 SAMPLE, PREVIEW, ANALYSIS, CATALOG, EXIT 

15 ? 

e. EXIT, Command; EX ) 

Description 

This command will simply end the program, with no intervention from the 

user. 

Notes 


If DAISE is mistakenly exited, type: 

BRUN DAISE. OBJ, Dl (cr) 
to start the program again. 

D. Second Level Subroutines 
Description 

SET-UP sets the analysis mode and loads the user specified data from 
disk into the APPLE® monitor to be ready for computation and analysis. User 
decisions are very important in this subroutine, therefore, the user is 
advised to read the computation algorithm section in this manual, preview 
the data, have a clear idea of what channel(s) are desired, what the inter- 
esting part of the data is, and how many averages are needed before 
attempting to enter this subroutine. 

User intervention 

1 First, input a ref-channel (subfile) name. This is done vfoether you wish 
to have an auto or cross-correlation analysis. 

2 Second, determine the option of AUTO-CORRELATION (AC) or CROSS- 
CORRELATION (CC). An AUTO-CORRELATION option is the route to (AFTER 
COMPUTATION) AMPL. and POWER SPECTRUM and AUTO-CORRELATION function 
of your ref-channel you entered in 1. 

While a CC option is the route to AMPL. spec of the ref-channel, x-power 
spec, between ref-channel and x-channel (that you will specify next and a 
CROSS-CORRELATION of the two with selection of the CC option. 

3 You will be asked to input a x-channel (subfile) name. 
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4 The 1 K bytes of data is divided into 8 blocks 128 bytes each. At this 

point you have to determine the block at which you want the computation 

to start. Your answer should be in the range of 0 to 7. 

5 Now you have to specify the number of blocks you want to be operated on 

by one pass through the ocmpute subroutine. Your answer should be in the 
range of 1 to 8. Note that this corresponds to the number of times you 
want your spectral results to be averaged. Also note that a second pass 
through the computer subroutine will double the number of averages (in 
5). However this number is limited to 8 average in any case, i.e., 

if your answer in 4^ is 2 
if your answer in J> is 3 
you can compute twice only. 


Please read the description of COMPUTE subroutine to help understand the 
procedure. 

Example 

1 TYPE COMMAND SE, CO, DI, EX, FOR: 

2 SET-UP, COMPUTE, DISPLAY, EXIT 

3 ? 

SE (cr) 

4 INPUT REF-CHANNEL 

5 ? 

TEST-C1 (cr) 

6 TYPE COMMAND AC, CC FOR: 

7 AUTO-CORRELATION OR CROSS-CORRELATION 

8 ? 

CC (cr) 

9 INPUT CROSS CHANNEL 

10 ? 

TEST-C5 (cr) 

11 STARTING BLOCK (0 - 7)? 

12 ? 

0 (cr) 

13 HOW MANY BLOCKS (1-8)? 

14 ? 

2 (cr) 

15 TYPE COMMAND SE, CO, DI, EX, FOR: 

16 SET-UP, COMPUTE, DISPLAY, EXIT 
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17 ? 


Notes 


In case of CC option, ref-channel should be the one that is expected to be 
delayed in time while the cross-channel is the advanced one. 

Warnings 

- Once you give the command SE all previous spectral and correlation 
results are lost and initialized to zero. 

If your answer in line 12 . is x 
and your answer in line 14 is y 

then x+y has to be < 8, otherwise DAISE will take you back to line 11 . 

If the number of passes through COMPUTE subroutine is z then: 

(z*y) +x < 8, otherwise COMPUTE will not execute but rather takes you 
back to line _1^, 2 and 2 . above. 

b. COMPUTE , command: CO 

Description 

Most of the computations are done in COMPUTE deploying the FFT routine. 
After execution, the real and imaginary components of ampl. and power 
spectrtm or ampl. and cross-power spectrun are computed. Ihese components 
will be the basis for the phase and magnitude displays as well as the cor- 
relation function in the third level subroutines. 

No results are obtained from COMPUTE unless you go through SET-UP once. 
Every pass through COMPUTE causes the mmber of data blocks specified in 
SET-UP to be operated on and the results are averaged with the previous ones 
obtained by an earlier pass. There is no user intervention in this sub- 
routine. 

Example 

Specified by SET-UP: Starting block is 1 

number of blocks is 3 

Then, as shown below, the first pass in COMPUTE operates on blocks 1, 2, and 
3 while the second operates on blocks 4, 5, and 6. 
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Note that a third pass will cause the error message: ERROR:DATA BLOCKS ARE 

EXCEEDED to be printed and control will be transferred to the beginning of 


the second level, since there is only 1 remaining block. 



+ 

start of first 
pass 


+ 

start of second 
pass 


t 

end of second 
pass 


Notes 

Due to the slow speed of APPLE* software in computing arithmetic func- 
tions, one pass through COMPUTE, operating on one single block, takes 70 
seconds. As a result, one pass operating on 8 blocks at once takes 9 
minutes and 20 seconds! 


c. DISPLAY, command: PI 

Description 

Assuming that COMPUTE has been executed and the spectral results are 
ready, DISPLAY provides viewing of these results on the video screen as well 
as on the SILENTYPE® printer. The functions that could be displayed are 
listed in the description of ANALYSIS subroutine earlier. 


DISPLAY takes DAISE into the third level of subroutine by grouping the 
displays of results into 3 categories as follows: 

1 - Amplitude spectrun as phase and magnitude. 

2 - Power or cross-power spectrum as phase and magnitude. 

3 - Auto-correlation or cross-correlation function. 

Description of these subroutines are provided in the subsequent sections . 
User intervention 

Only one decision has to be taken at this stage, that is choosing one 
out of 3 subroutines: AMPLITUDE SPECTRUM, POWER SPECTRUM, or CORRELATION. 

These subroutines are executed in any order, and as many times as desired. 
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Example 


1 TYPE COMMAND SE, CO, DI, EX, FOR: 

2 SET-UP, COMPUTE, DISPLAY, EXIT 

3 ? 

DI (cr) 

4 TYPE COMMAND AS, PS, CR, EX, FOR: 

5 AMPLITUDE SPECTRUM, POWER SPECTRUM, CORRELATION, EXIT 

6 ? 

PS (cr) 

d. EXIT, command: EX 

Description 

This command exits you from the second level subroutines back to the 
first level. 

Example 

1 TYPE COMMAND SE, CO, DI, EX, FOR: 

2 SET-UP, COMPUTE, DISPLAY, EXIT 

3 ? 

EX (cr) 

4 TYPE COMMAND SA, PR, AN, CA, EX, FOR 

5 SAMPLE, PREVIEW, ANALYSIS, CATALOG, EXIT 

6 ? 

E. Third Level Subroutines 

a. AMPLITUDE SPECTRUM, command: AS 

Description 

AMPLITUDE SPECTRUM generates the phase and magnitude of the amplitude 
spectrum of the reference channel specified in SET-UP. The viewing options 
are: 


1 - Phase display 

2 - Magnitude display 

3 - Printing real and imaginary values of spectral components on the 

video screen. 
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These 3 options are, in turn, 3 subroutines that are called with commands as 
the previous ones. The order in which these routines are called is irrele- 
vant and is left to the user 1 s convenience. 

1 - PHASE 

Command: PH 

In this case the first 128 components of the phase are drawn with an argu- 
ment of + 180' to -180'. Underneath the display, a line is printed contain- 
ing the name of the reference-channel. Next the user is asked whether he 
wishes to print the display using the SILENTYPE® printer or not. His answer 
must be Y (yes) or N (no) . If it is a Y, then the SILENTYPE® prints a 
heading followed by the display, and the control is transferred back to the 
3 options above. If the answer is N, then control will be transferred 
directly to level 3. 

Example 

1 TYPE COMMAND AS, PS, CR, EX, FOR: 

2 AMPLITUDE SPECTRUM, POWER SPECTRUM, CORRELATION, EXIT 

3 ? 

AS (cr) 

4 TYPE COMMAND PH, MG, R & I, EX, FOR: 

5 PH-DSPLY, MG-DSPLY, REAL & IM. , EXIT 

6 ? 

PH (cr) 

(NOW THE PHASE IS PLOTTED) 

7 XXXXPHASE PLOT: TEST-C2 XXXX 

8 WOULD YOU LIKE TO PRINT THE DISPLAY (Y/N) 

9 ? 

Y (cr) 

THE SILENTYPE® WILL PRINT: 

PHASE PLOTXXX TEST-C2XXX 4 AVE 
AND THEN THE DISPLAY IS PRINTED. 

10 TYPE COMMAND PH, MG, R & I, EX, FOR 

11 PH-DSPLY, MG-DSPLY, REAL & IM., EXIT 
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12 ? 


Notes 


An invalid command in line 6 or 9 will be refused and the question will be 
asked again. 

Example 

1 TYPE COMMAND PH, MG, R A I, EX, FOR 

2 PH-DSPLY, MG-DSPLY, REAL A IM. , EXIT 

3 ? 

MG (cr) 

4 SUPPRESS SPECTRAL COMPONENTS (Y/N)? 

5 ? 

Y (cr) 

6 WHICH COMPONENT (0-127)7 

7 ? 

8 SUPPRESS SPECTRAL COMPONENTS (Y/N)? 

9 ? 

N (cr) 

10 FULL, HALF, QUARTER DISPLAY (F/H/Q)? 

11 ? 

H (cr) 

12 SCALING FACTOR? 

13 ? 

1 (cr) 

14 SCALING FACTOR? 

15 ? 

10 (cr) 

(NOW THE MAGNITUDE IS PLOTTED) 

16 XXX MG. PLOT: TEST-C3 XXX 

17 WOULD YOU LIKE TO PRINT THE DISPLAY (Y/N)? 

18 ? 

N (cr) 

19 TYPE COMMAND PH, MG, R A I, EX, FOR: 

20 PH-DSPLY, MG-DSPLY, REAL A IM., EXIT 
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21 ? 


2 " magnitude 

Command: MG 


This option plots for you the first 128 components of the magnitude plot-. 
The vertical axis represents a relative scale while the horizontal one is 
frequency. 


Before plotting the magnitude, the user has to answer the following 
questions: 

i - suppress spectral components (Y/N) ? 

This allows suppression of relatively large components 
(i.e. DC comp.) in order to have a better view of the rest of 
the spectrum by changing the scaling factor. If your answer is 
yes, then the user will be asked to specify the number of the 
components which can only be suppressed, one at a time. 

ii - Full, half or a quarter display (F/H/Q)? 

This allows changing the increment value on the frequency 
axis by displaying 128, 64 or 32 cong>onents. 

F corresponds to 128 components 
H corresponds to the first 64 components 
Q corresponds to the first 32 components 

iii - Scaling factor? 

This input could be any real number by which all components are 
divided before they are plotted. 

If the scale factor is too small the question will be repeated 
until the scaling factor is large enough for the components to 
fit on the APPLE* screen. 


After answering these questions the magnitude is plotted. As in the PHASE 
plot information is printed underneath followed by the option of printing 
the display. 

REAL AND IMAGINARY 

Command: R & I 

This option allows printing the real and imaginary value of any spectral 


I 
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component on the video screen. This command prints the values of one coo* 
ponent at a time. 

Example 


1 TYPE COMMAND PH, MG, R & I, EX, FOR: 

2 PH-DSPLY, MG-DSPLY, REAL & IM. , EXIT 

3 ? 

R A I (cr) 

4 WHICH COMPONENT (0-127)? 

5 ? 

2 'cr) 

6 REAL PART: 2356.712 

7 IM.PART: -175.689 

8 TYPE COMMAND PH, MG, R & I, EX, FOR: 

9 PH-DSPLY, MG-DSPLY, REAL & IM. , EXIT 

10 ? 

Notes 

- The frequency resolution on the horizontal axis, in both PHASE and MAGNI- 
TUDE, is determined by the sampling frequency specified in the SET-UP. 

Th« relation i. A £ - «-P lin » *.”■ 

256 

- Along with PH, MG and R A I you have the EX command. This command takes 
DAISE back to level 3 subroutines. 

- Suppressing spectral components in MAGNITUDE means setting their values 
to zero only on the display so that in subsequent passes through MAGNI- 
TUDE these components will still be suppressed while a second call of 
AMPLITUDE SPECTRUM will regenerate them. 

b. POWER SPECTRUM, command: PS 

Description 

POWER SPECTRUM generates the phase and magnitude of the power spectrum 
of the reference-channel if the auto-correlation option was specified in the 
SET-UP, or it generates the phase and magnitude of the cross-power spectrum 
of the reference-channel and cross-channel if the cross-correlation option 
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was specified. 


Next, it leaves the user with the sane options listed in the AMPLITUDE 
SPECTRUM description. From that point on, everything mentioned there is 
true in this case also and need not be repeated here. 

Notes 

The comments printed underneath the video displays and above the print- 
ed display* will refer in this case, to the power or cross-power spectrum 
rather than the AMPLITUDE SPECTRUM. 

c. CORRELATION, command: CR 

Description 

This subroutine uses Che results of the power, or cross-power, spectrum 
to generate the auto-correlation, or cross-correlation function by one pass 
through the inverse FFT. 

After generation, the correlation function is displayed and could be 
printed by the printer. 

The vertical axis represents a relative sca?e of the function while the 
horizontal one represencs the delay time t. 

The number of points plotted is 256 with a time resolution determined 
by the sapling frequency as: 

At - 

r s 

User intervention 

1 At the very beginning, the user has the option of suppressing spectral 
components in the power spectrum before generating the correlation- func- 
tion. 

This is included here to prevent undesirable components (i.e. noise com- 
ponents) from effecting the correlation results. 

2 Once the function is generated, the axis will be drawn and the user is 

asked to input a scaling factor for the display. 

3 A comment is printed underneath the display specifying the kind of func- 

tion and the involved channel(s). Next the question: 


' 22 


Diaplfy CORR-FN again (Y/N)? 


Is aaked. This allows viewing Che function on a different scale, and 
could be repeated ae many times as desired. 

4 Then the option of printing the display on the printer is offered. 

Example 

1 TYPE COMMAND AS, PS, CR, EX, FCR: 

2 AMPLITUDE SPECTRUM, POWER SPECTRUM, CORRELATION, EXIT 

3 ? 

CR (cr) 

4 SUPPRESS ANY SPECTRAL COMPONENTS (Y/N)? 

5 ? 

Y (cr) 

6 WHICH COMPONENT (1-238)? 

7 ? 

3 (cr) 

8 SUPPRESS ANY SPECTRAL COMPONENTS (Y/N)? 

9 ? 

N (cr) 

(AFTER 70 SECONDS OF COMPUTATION THE AXIS ARE PRAWN) 

10 SCALING FACTOR? 

11 ? 

8 (cr) 

(THE FUNCTION IS PLOTTED) 

12 XXX CROSS -CORR: TEST-C2 CROSS TEST-C7 XXX 

13 DSPLY CORR-FN AGAIN (Y/N)? 

14 ? 

Y (cr) 

15 SCALING FACTOR? 

16 ? 

4 (cr) 

(THE OLD FUNCTION IS DELETED AND THE NEW ONE IS PLOTTED) 

17 DSPLY CORR-FN AGAIN (Y/N)? 

18 ? 

N (cr) 


% 
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19 WOULD YOU LIKE TO PRINT THE DISPLAY (Y/N)? 

20 ? 

N (cr) 

21 TYPE COMMAND AS, PS, CR, EX, FOR: 

22 AMPLITUDE SPECTRUM, POWER SPECTRUM, CORRELATION, EXIT 

23 ? 

Notes 

“ In this subroutine, suppressing spectral components means setting them to 
zero everywhere in the program and they cannot be generated again. 

i.e., A suppressed component in one pass through CORRELATION will still 
be suppressed in any subsequent passes, as well as in any power spectrum 
display. 

Also note that the DC component is suppressed automatically before computing 
the correlation-function. 

- Once the user answers "N" to the question in line 17 of the example 

above, the CORRELATION subroutine is exited. Therefore, a second view of 
the correlation function will require a second pass through the sub- 
routine which in turn means waiting 70 seconds. 

d. EXIT, command: EX 

Description 

This command takes DAISE back to level 2 subroutines. 

F. Data Management 

Data is stored and transferred between the APPLE* monitor and disk 
drive 2 at three major subroutines: SAMPLE, PREVIEW and SET-UP. 

a. SAMPLE Data 

In this case data is sampled at the A/D converter and stored directly 
in the APPLE* RAM as 8 blocks (channels), 1024 bytes each, stored at the 
following locations: 

Channel 1 starts at 7600 HEX (30208 DEC) 
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Channel 2 starts at 7AOO HEX (31232 DEC) 

Channel 3 starts at 7E00 HEX (3256 DEC) 

Channel 4 starts at 8200 HEX (33280 DEC) 

Channel 5 starts at 8600 HEX (34304 DEC) 

Channel 6 starts at 8A00 HEX (34328 DEC) 

Channel 7 starts at 8E00 HEX (36352 DEC) 

Channel 8 starts at 9200 HEX (37376 DEC) 

These blocks are transferred to disk drive 2 as 8 binary subfiles and the 
memory space (7600 to 9600 HEX) is freed. 

b. PREVIEW Data 

In the PREVIEW subroutine one subfile (1 k byte) is loaded from disk 
drive 2 into the RAM and occupies: 

8E00 to 9 IFF HEX (36352 to 37375 DEC) 

Once PREVIEW is exited this 1 K block of data is not accessed by any means. 

c. SET-UP Data 

In this case one subfile (reference-channel) is loaded from the 
diskette and stored in locations: 8E00 to 91FF EX(36352 to 37375 DF.C) and 

in the case of cross-correlation the x-channel is stored in 9200 to 95FF HEX 
(3 7376 to 38399 DEC) 

These 2 blocks of data are accessed by COMPUTE subroutine only. 

Notes 


- Memory locations 7600 to 95FF are used only for data management through- 
out DAISE. 

- All spectral and correlation results are not stored permanently anywhere 
and will be lost unless the displays are printed, 

III. COMPUTATION ALGORITHM 

The algorithm used in DAISE for the computation of AMPLITUDE SPECTRUM 
and POWER SPECTRUM is based on two algorithms as described by Cooley et al. 
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[2] and Rader [ 3 J . The reader is advised to review the papers for clear 
under st and lug. For our purposes a summary of the two algorithms is 
presented. 

A. First Reference 


From reference [2j we can compute the Fourier transform of 2 sets of 
real data in one pass through a DFT subroutine as follows: 

Using the linearity property we see that if x(n) and y(n) are real sequences 
such that 

X(K) is the transform of x(n) and 
Y(K) is the transform of y(n) 
and if we form s(n) =* y(n) + i x(n) 

Then $(n) has the transform 

S(K) - Y(K) + i X(K) (1) 

and S*(N-K)“ Y(K) - i X(K) (2) 0 < K < j 


where * designates complex conjugate and N is the total number of points. 
Now solving (1) and (2) for Y(K) and X(K) we get: 

Y(K) - 1/2 [S*(N-K) + S (K)] 

X(K) - 1/2 [ S*(N-K) - S(K)] 

Thus, Y(K) and X(K) can be computed by one pass through the DFT with a few 
extra additions and subtractions. 


B. Second Reference 

With respect to perform correlation operations as described by Rader 
[3], let Y(N) be an N points sequence and construct X(N) such that: 


x(n) 


y(n) for 0 < n < — 


0 for — < n < N n£ [0,N) 
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ORIGIN At HAU* • . 

If we compute the DFTs of Y(N) and X(N): OF POOR QUALITY 

Y(K) - DFT { y(n)} 

X(K) - DFT { x(n)} 

and for the product W(K) ■ Y*(K) X(K) we have the DFT, W(K) of the sequence 
N/2 -1 

w(m) ■ £ x(n)*y(n + m) (2) 

n“0 

where m is a lag index. 


It is easily seen that w(m) is nothing but the auto-correlation function of 
y(n) that, having W(K) computed, can be generated by on inverse DFT. 

IF y(n) was a continuous series the above procedure can be repeated as many 
times as desired, considering N points at a time and averaging. 


Note that Rader paper is more involved and does not stop at this point, but 
the above discussion is enough for our need. 

C. Algorithm Description 

Our algorithm is the merit of both papers. Let y(n) and x(n) be 2 
real sequences of length N. The second ^ points of x(n) are padded with 
zeroes as follows: 


1) 


2 ) 


If auto-correlation is desired, 

|y(n) for 0 < n < n/2 


x(n) ■ 


| 0 for N/2 < n < N 
If cross-correlation is desired, 

I x(n) for 0 < n < N/2 


x(n) = 


0 for N/2 < n < N 
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then s(n) ■ y(n) + j x(n) is formed and its Fourier transform S(K) * a(K) + 
i b(K) is obtained by one pass through the FFT routine. 

a. Amplitude spectrum computation 

We have seen in the first procedure that 

Y(K) * 1/2 [ S*(N-K) + S (K) ] 

or using the results of the FFT routine, 

Real part of Y(K) - 1/2 [ a(N-K) + a(K)j » Real part of Y(N-K) 

and imaginary part of Y(K) ■ 1/2 [ -b (N-K) + b(K) ] = -Imaginary part of 
Y(N-K) 

The last two expressions are the ones used to evaluate the AMPLITUDE 
SPECTRUM. 

b. Power or Cross-Power Spectrum Computation 
We saw that Y(K) =* 1/2 [ S*(N-K) + S(K)] 

and its complex conjugate would be 
Y*(K) * 1/2 [ S*(N-K) + S*(K) ] 

so that the POWER SPECTRUM W(K), in the Rader paper can be written as: 

W(K) = Y*(K) X(K) - 

1/2 [ S(N-K) + S*(K)] • 1/2 [S*(N-K) - S (K) ] . 

Evaluation of the above expression in terms of the output of the FFT 
routine, S(K) * a(K) + i b(K) yields, 

Real part of W(K) * Real part of W(N-K) * 

1/2 [ a(K) b(N-K) + a(N-K) b(K)] 

and imaginary part of W(K) e - Imaginary part of W(N-K) 

- 1/4 [a 2 (N-K) - a 2 (K) + b 2 (N-K) - b 2 (K) ] 

where 0 < K < ^ - 1 

Again the last two expressions are the ones used for evaluation. 

c. A lgorithm Summary 

Step 1 Load data sequences y(n) and x(n), and then zero pad the second 
W/2 points of x(n) . 
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Step 2 Form s(n) as y(n) + i x(n) and pass it to the FFT routine to 
get 

S(K) - a(K) + i b(K) 

Step 3 Compute the AMPLITUDE SPECTRUM of y(n) and the power (or cross- 
power) spectrum as specified earlier. 

Step 4 Repeat the above steps as many times as desired and average the 
results in Step 3. (In DAISE the niznber of averges is limited 
to 8) . 

Step 5 Pass the power (or cross power) spectrum results to the inverse 
FFT to obtain the auto (or cross) correlation function. 


Notes for the User 


1 - In DAISE y(n) has been called the reference-channel while x(n) is the 
cross-channel. 


2 - Amplitude spectrum of y(n) is computed whether the user specifies the 
auto-correlation or the cross-correlation option. 


3 - Our sequences (subfiles) are 1024 points long, divided into 8 blocks, 

128 points each. The value of N is 256 points and 

y - 128 - 1 block. 

4 - DAISE deploys a 256 point FFT routine so that by one pass through the 

first 3 steps of the algorithm two blocks of y(n) and 1 block of x(n) 
(second block is padded with zeros) are operated on. This leaves us 
with 8 possible passes or averages. 


l .e. 


First 

average : 

y(n) * 

Second 

average : 

x(n) - 
y(n) = 

and so 

on. 

x(n) ■ 


first and second block} 
first block and 128 zeros} 
second and 3rd block} 
(second block and 128 zeros} 


5 - In SET-UP subroutine, when the user answers the question "HOW MANY 

BLOCKS (1-8)?" He is specifying the number of averages desired in one 
pass through COMPLETE subroutine. 


ORIGINAL PAtic 
OF POOR QUALITY 
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APPENDICES 
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ORIGINAL PALL E!s 

OF POOR QUALITY 


A. Memory Map 


0000 (80000) 
2048 (80800) 

6063 (817AF) 
8192 (82000) 

16364 (84000) 

26666 (8682A) 
30208 (87600) 

38400 (89600) 

45056 (8B000) 
49151 (8BFFF) 



Note 

In the above configuration DOS, the disk operating system has been moved 
to the extra 16 K memory card. 


B. Time-Frequency Relations 


F g ■ sampling frequency 

Af - frequency resolution (spectral plots) 


At ■ time resolution (correlation plots) 


T ■ total sampling time (N points) 


relations: 


At 


Af 


s 

F 

_t 

N 


r - — 


Af 


Examples: a) F » 128 Hz, N » 256 


At - 


Af 


1 

128 

128 

256 


.0078 seconds 


.5 Hz 


r ■ “J * 2 seconds 


b) F - 51.2 Hz, N - 256 


At 


Af 


51.2 

51.2 

256 

1 


.02 seconds 
.2 Hz 
5 seconds 
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C. DAISE Error Messages 

a. "Out of Range” Errors 

- Input in line 1080 

Message: "ERROR: QUARTER RANGE IS 1 TO 4, INPUT AGAIN" 

- Input in line 1870 

Message: "ERROR: YOUR INPUT RANGE IS 0 TO 7, INPUT AGAIN" 

- Input to line 1890 

Message: "ERROR: YOUR INPUT RANGE IS 1 TO 8, INPUT AGAIN" 

- If input in line 1870 plus input in line 1890 is greater than 8 

Message: "ERROR: DATA RANGE IS EXCEEDED!" 

- Input in line 5513 or line 5724 

Message: "ERROR: COMPONENT RANGE IS 0 TO 12, INPUT AGAIN" 

- Input in line 5885 

Message: "ERROR: COMPONENT RANGE IS 0 TO 128, INPUT AGAIN" 

b. DOS Related Errors 

If a DOS error occurs while DAISE is storing data on disk after sampl- 
ing or loading data from disk for PREVIEW or ANALYSIS the corresponding 
error message will be printed and control will be transferred to the first 
level subroutines. 

Error Messages: 

- "ERROR: DISKETTE IS WRITE PROTECTED" 

- "ERROR: FILE NOT FOUND" 

- "ERROR: DISKETTE VOLUME MISMATCH" 

- "ERROR: I/O ERROR" 

- "ERROR: DISK IS FULL" 

- "ERROR: FILE IS LOCKED" 

- "ERROR: BAD FILE NAME" 

- "ERROR: UNKNOWN", (if error is not one of the above) 
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D. Flow Chart 


Start 


ORIGINAL PAul! 1$ 
OF POOR QUALITY 



To second level subroutines 

next page 
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From spectrum 
subroutines 


IT 














E. Complete Listing 


ccv.-:**.: • • t 

OF POOR Q’JAUfY 

iLIST 

16 REM 
30 DS * 

36 ONERR SOTO 400 

100 DIM M< 127 )s DIM K 12?) 

103 DIM ZR< 128 ): DIM 2K128) 

105 DIM Vft< 123): DIM VI<128> 

187 DIM AR<255): DIM AI<255> 

120 PRINT D*“BLQAO ADC. 0,0*300,01" 

140 PRINT "TVPE CCPlMAND SA, PR, AN, CA, EX, FOR: " 

150 PRINT “ SAMPLE, PREOIEW,ANALYSIS»CATALOG,EX‘IT" 

1S0 I NPUT A* 

1?0 IF A* = “SA" THEN GOSUB 500 
180 IF A* * -PR“ THEN 60SUB 1000 
190 IF A* * "AN" THEN GOSUB 1500 

209 IF A* * -CA" THEN TEXT : PRINT DS"CATAL06,D2": SOTO 140 

210 IF A* a "EX" THEN SOTO 238 
220 SOTO 140 

230 END 

480 REM DOS ERROR DETECTION 
485 DED a peek <222) 

410 ON DED SOTO 415,415,415,420,415,425,430,435,440,445,450,415,415,415,4 
15 

415 PRINT “ERROR: UNKNOWN": SOTO 148 

420 PRINT "ERROR: DISKETTE IS WRITE PROTECTED": SOTO 140 

425 PRINT “ERROR: FILE NOT FOUND": SOTO 140 

438 PRINT "ERROR: DISKETTE OOLUME MISMATCH": SOTO 140 

435 PRINT “ERROR: I/O ERROR": SOTO 140 

448 PRINT "ERROR: DISK IS FULL"? SOTO 140 

445 PRINT “ERROR: FILE IS LOCKED": SOTO 140 

45© PRINT "ERROR: BAD FILE NAME": SOTO 140 

508 REM SAMPLE 

505 PRINT “INPUT SAMPLING FREQUENCY CONSTANT" 

518 INPUT SFD 

515 IF SFD < 1 THEN PRINT "ERROR: CONSTANTCAN NOT BE < 1, INPUT AGAIN":. GOTO 
510 

520 POKE 737, SFD 
530 CALL 768 

535 PRINT "TYPE FILE NAME: “ 

548 INPUT FI* 

550 Cl* * FI* + “-Cl" 

551 C2* = FI* + "-C2" 

552 CSS = FIS + "-C3" 

553 C4S = FT.* + "**C4" 

554 CSS * PI* + “-C5" 

555 CSS a FIS + "-C6" 

55S C7S » FI* + U -C?“ 

557 CBS = FIS 4 "-C8" 

560 PRINT OS ; " BSAUE" ;C 1 S ; " , AS7600 , LI 024 , D2 " 

570 PRINT OS;"B$AUE" ;C2S;",A$?A00,L1024,D2" 

588 PRINT D*;"BSAUE" ;C3S;“ ,A*?E00, LI 824,02" 

590 PRINT D*5“BSAUE" ;C4*;“ ,A*8200,L1024,D2 ,, ‘ 


39 


600 

616 

620 

630 

640 

700 

705 

710 

715 

720 

300 

310 

820 

S30 

335 

340 

345 

850 

355 

857 

370 

890 

900 

910 

915 

928 

925 

930 

935 

940 

945 

950 

955 

960 

965 

978 

971 

972 

973 


PRINT OS*"BSAUE";C5*i"*Pf8S80*L1024*O2“ 

PRINT W;“BSOUE-iC6«i“,0*8fi00,L1024,O2 u 
PRINT D$*"BSAUE"*C?$;"*A$8E00*L1U24,O2“ 

PRINT OS;-BSOUE-iC8S;“,R*9200,U024,D2 , ‘ 

RETURN 

PRINT “WOULD YOU LIKE TO PRINT THE DISPLAY* Y/N)?- 
INPUT PD$ 

IF PQS * “V" THEN 6GSUB 300: RETURN 
IF PD* = "N" THEN RETURN 
SOTO 700 

REH PRINT DSPLY 
PR# 1 

L* * REH CTRL Q 
PRINT 

POKE - 12529*255 
PRINT 

POKE - 12524*0 
PRINT 

POKE - 12528*7 
Tfl * T / 4 


S-Zi 

0^ POOR QUALITY 


IF K* * “CR“ THEN GOTO 955 
IF KS * “AS* THEN SOTO 940 
IF E * 99 THEN SOTO 925 
IF F* * “PH“ THEN GOTO 920 

PRINT “X-POHER***“RCHS" CROSS "XCH*"» U TR" PU63H»SCPLE*“fl: SOTO 970 
PRINT “X-POHSfc»*“RCH*" CROSS "XCH*“***-TA“ AUS“: SOTO 97 9 
IF FS * “PH“ THEN SOTO 935 

PRINT “POWER SPEC.***-RCH*"***“TA" AUS***SCALE=“A: SOTO 979 
PRINT "PHASE OF POWER SPEC«hs-“RCH*“*-**“TA m PI'S": SOTO 970 
IF Ff = “PH“ THEN SOTO 950 

PRINT "PHPL. SPEC.*-**“RCH*“***“TA“ AUG***SCAL£a"A: GOTO 976 
PRINT "PHASE OF AMPL. SPEC*^“RCH*"***"TA" PUS": GOTO 978 
IF E = 99 THEN GOTO 965 

PRINT "X-CORR***"RCH*" CROSS “XCHS"*~^"TA“ AU6***5CALE*“U: SOTO 970 

PRINT "CORR***"fiCH*-***“TR‘‘ AU6***SCALE*‘ , U 

PRINT 

PRINT 

PRINT 

PRINT L* 


974 

PR# 0 



975 

RETURN 



1000 

REH 

PREUIEW 


1050 

PRINT 

"DESIRED FILE PND 

CHANNEL?" 

1055 

INPUT 

RCHS 


1070 

PRINT 

"DESIRED QUARTER* 1-4 >?" 

1080 

INPUT 

Q 


1082 

IF Q -■ 

: 1 OR Q > 4 then 

PRINT "ERROR: QUARTER RANGE IS 


PSPIN’ 

1 : SOTO 1080 


1085 

SOSUB 

4000 


1030 

GOSl© 

4200 


1035 

M = 0 



1100 

GOSL© 

M 

8 

© 


1110 

PRINT 

"****• REAL DATA: 

"RCHf" , QR- M Q" ***»" 


1 160 RETURN 

1200 REH DSPLY REAL OPTO OR CORR-fN 

1220 H6F: : HCOLOR* 7 

1225 HPLOT 0#0 TO 0*159 

1230 FOR V = 159 TO © STEP - 32 

1235 HPLOT 0*Y TO 3*Y: NEXT 

1240 IF W * 0 THEN 60 TO 1265 

1245 HPLOT 0,79 TO 255*79 

1250 FOR X » 0 TO 255 STEP 40 


1 TO 4, INPUT 
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1255 HPLGT X.73 TO X.76 
1260 NEXT X: 60T0 1302 
1265 HPLGT 0.159 TO 255.159 
1270 FOP. X * 6 TO 255 STEP 46 
1275 HPLOT X.153 TO X.156 
1260 NEXT X 

1302 PRINT “SCALING FACTOR?" 

1304 INPUT 1.1 

1310 FOR X - 0 TO 255 

1320 V * 359 - H - ARC X > * .623 / U 

1322 IF V > 159 THEN GOTO 1302 

1324 IF V < 9 THEN SOTO 1302 

1330 HPLOT X + l.Y 

1346 NEXT 

1 DCTI ©W 

1500 REM SPECTRUM ANALYSIS 
1550 PRINT "TYPE COMMAND SE, CO. DI. EX. FOR: 
1560 PRINT “SET-UP. COMPUTE. DISPLAY, EXIT-* 
157© INPUT Of 

1580 IF C$ = “SE" THEN SOSUB 1650 

1590 IF CS = “CO- THEN 60SUB 2166 

1S00 IF CS * ”01“ THEN 60SUB 5800 

1610 IF Cl = “EX- THEN SOTO 1636 4 

1620 SOTO 1550 

1630 RETURN 

1650 REH SET-UP 

1708 FOR X * 0 TO 128 

1718 2R( X > * 0 

1720 ZKX > * 0 

1730 YR( X > * © 


L. « . < l . . ; 'Li it 1 j 

or poor quality 


ii 


1740 YI(X > * 0 
1750 NEXT 

1760 PRINT "INPUT REFRENCE CHANNEL- 
1770 INPUT RCHS 

1730 PRINT “TYPE COMMAND AC, CC, FOR:" 

1790 PRINT “AUTO-CORRELATION, CROSS-CORRELATION" 
I860 INPUT BS 

1810 IF BS = “AC“ THEN E * 99: SOTO I860 
1820 IF BS = “CC" THEN SOTO 1840 
1830 SOTO 1730 
1840 PRINT "INPUT X-CHANNEL" 

1850 INPUT XCHS 


I860 PRINT "START I NS BLOCK* 0-? >?“ 

1870 INPUT S 

1875 IF S < 0 OR S > 7 THEN PRINT “ERROR: YOUR INPUT RANGE IS 0 TO 7. IMP 
UT AGAIN": GOTO 1860 
1 830 PRINT “HOW MANY BLOCKS* 1-8 >?“ 

1890 INPUT h 

1892 IF H < 1 OR H > 8 THEN PRINT “ERROR: YOUR INPUT RANGE IS 1 TO 8, IN 
PUT AGAIN": GOTO 1880 
1895 T * 2 * H 

1900 IF 03 + H > > = 9 THEN PRINT “ERROR: DATA RANGE IS EXEEDED ‘ " : GOTO 
1 860 

1910 GO SUB 4000 

1920 IF E - 93 THEN GOTO 1948 

1930 GOTO 1960 

1940 RETURN 

1950 REM LOAD CROSS CHANNEL 

1 96.0 PR I NT OS ,* " BLO AD “ i XCHS ; " , A$9£00 , D2 " 

1970 SOTO 1940 
2100 REH COMPUTE 

2150 P = 256 
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2160 IF <6 + H> > * 9 THEN PRINT “ERROR: 

550 

2178 J = 36352 + < S - 1 ) * 123 

2180 FOR V e l TO H 

2130 J = J + 123 

2200 FOR I ~ 0 TO 255 

2210 AR<I) = PEEK CJ + I) 

2230 NEXT 

2240 IF E - 33 THEN K * J: SOTO 2260 
2256 K = J + 1024 
2266 FOR 1=0 TO 12? 

2276 AI< I > * PEEK'. < K * I ) 

2286 fll< I + 128 > * 0 
2285 NEXT 
2280 FRO * 1 
2380 6GSUB 3800 

2382 2R<0> = ZRC0) + «FK0> * AI<0> 

2304 VR<0) * VR(0> -+ fiR<0> 

2310 FOR X * 1 TO 128 

2320 R = < PFK X ) > * (flICP - X>> 

2330 S = <GR(P - X)> * <AICX)> 

2340 2R( X > * 2R< X > + < R + S > ✓ 2 
2350 R = <AR<P - X)> ^ 2 - <OROO> ^ 2 
2360 S = (flKP - X)> - 2 - <flI<X>> - 2 
2370 ZI<X> = 2KX ) - <R ♦ S> ✓ 4 
2380 R = AR< P - X> + AR<X> 

2390 S = fil(X) - AI<P - X) 

2480 VRvX) = VR<X> + R / 2 

2410 VKX) s yi(X) + S / 2 

2420 NEXT X 

2430 NEXT V 

2435 T * 2 * T 

2440 S “ S + H 

2450 RETURN 

3000 REM FFT ROUTINE 

3100 M a 8 

3110 N * 2 ^ H 

3120 HI = 0:N2 = N - 1 

3130 FOR N3 * 1 TO N2 

3140 M4 — N 
3156 N4 * N4 / 2 

3160 IF N1 + N4 > N2 T HEN SOTO 3150 
317© Ml = N1 ~ I NT HI / N4) * N4 + N4 
3180 IF N1 < N3 THEN SOTO 3256 
3190 Tl = Afi<N3> 

3200 AR< N3 > = flRi.Nl) 

3210 AR< HI > = Tl 
3220 T2 = AK.N3) 

3230 AKN3) = flKNl) 

3240 flKNl) = T2 
3250 NEXT N3 

3260 REH DO THE COUPLE*: TRANSFORM 

3270 N4 = 1 

3280 N6 - 2 * N4 

3290 FOP N3 = 0 TO N4 - 1 

3300 fl = FP.D * N3 * 3.141592? / N4 

3310 C = COS «: fi » 

3320 S = SIN (A) 

3330 FOR N7 = N3 TO N - 1 STEP N6 

3340 N8 = N? + N4 

3350 Tl * C * AR<N8> - S * AI(N8> 

3360 T2 * C * AICN8) + S * AR<N3). 


DATA RANG IS EXEEDED!": 


ORtGm 'AL 

°p poor 


Quality 


WTO 1 
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3370 ARCN8) * ARC N7 ) - Tl 
3330 AICN8) * AICN7 ) - T2 
3330 ARC N? ) * ARC N? > + Tl 
3400 file N7> = AICN? > + T2 
3416 NEXT N7 
3420 NEXT N3 
3430 N4 » N6 

3440 IF N4 < N THEN SOTO 3230 
3442 FOR X * 8 TO 255 

3444 ftft< X > = ARCX) ✓ 32 

3445 file X > * AICX) s 32 
3448 NEXT X 

3450 RETURN 

4000 REM LOAD REF. CHANNEL 

4650 PRINT t»i"BLOAD*‘;RCHS;“,AS8E00#D2 M 

4080 RETURN 

4200 REN LOAD REAL ARRAV WITH ONE QUARTER 
4245 X = 36352 

4250 ON Q SOTO 42C3,42?0,4288,4230 

4260 J = X: SOTO 4300 

4270 J = X + 256: SOTO 4300 


w; 


OF POOR QUALliV 


4280 -J * X + 512: GOTO 4300 
4290 J * X + 768 
4380 FOR I - 0 TO 255 
4310 AR< I > = PEEK CJ + I) 

4328 NEXT I 
4330 RETURN 

4745 IF SOt = “Q“ THEN SO = 8: SOTO 5748 
5000 REH DISPLAY 

5050 PRINT “TYPE COWHAND PS* AS, CR, EX, FOR: “ 

5860 PRINT -POWER SPEC# AHPL. SPEC, CORRELATION, EXIT- 
5070 INPUT KS 


5080 IF K$ = "PS* THEN SOTO 5155 

5090 IF KS * "AS" THEN SOTO 5260 

5100 IF KS * '*CR“ THEN GOTO 5850 

5110 IF KS * “EX* THEN GOTO 5140 

5120 SOTO 5050 
5140 RETURN 

5150 REM POWER SPECTRUM 
5155 REM 


5160 FOR X * 0 TO 127 

5165 MCX) * SQR <<ZRCX» 2 * CZKX>> * 2) 

5170 IF ABS CZICX>> < T THEN SOTO 51S8 
5175 SOTO 5190 

5188 IF ZRC X ) < = - 7 THEN PCX) * 3.1415927: SOTO 5236 

5185 PCX) = 9: SOTO 5230 

5130 IF ABS < ZRC X ;• > < T THEN SOTO 5200 

5195 SOTO 5210 

5200 IF ZICX> > = T THEN P<X> • .5 * 3.1415327: GOTO 3230 


5205 PCX) = - .5 * 3.1415327: SOTO 5230 
5216 IF ZRC X > > » T THEN PCX) = ATN CZICX) 


/ < ZR< X > 


. 000001 > 


230 


5215 IF ZleX) > = T THEN PCX,' = 3.141592 - 

5230 

5220 PCX) = - 3.1415327 - ATN CZICX) •••" < - 


5230 

NEXT X 


5240 

SOSUB 5400 


5250 

GOTO 5050 


52S0 

REH REF-CH 

SPECTRUM 

5265 

REH 


5267 

FOR X * Q TO 

127 

5270 

MCX) = SQR C 

< yp< x > ■ 


ATN e ZIC X ) 

ZRC X > > ) 


7QC 


SOTO 5 
!>':■>: SOTQ 
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5275 IF BBS CVICX)) < T THEN GOTO 5285 ORi^h'^i r , . . p . 

5280 GOTO 5295 OF Pq Gv r .; . 

5285 IF VROO < * - T THEN POO = 3.1415927: GOTO 5349 ’ v ‘ 4, *Y 

5290 PC X > * 0: GOTO 5340 

5295 IF BBS CVP (OO) < T THEN (GOTO 5305 

5300 GOTO 5315 

5305 IF VIOO > * T THEN PCX) * .5 * 3.1415927: GOTO 53*8 
5310 WX) = - .5 * 3.1415927: GOTO 5340 

5315 IF VROO > * T THEN PCX) = BTN CVICX) ✓ CVRCX) + .000001 >>: GOTO 5 
340 

5320 IF VI(X) > a T THEN PCX) = 3.1415927 - BTN CVICX) ✓ < - VRC X)>>: GOTO 
5340 

5325 PCX) = - 3.141592 - BTN CVICX) ✓ C - VRC X)>): GOTO 5340 

5340 NEXT 

5350 GOSUB 5400 

5366 GOTO 5050 • 

5400 REM PHASE & MAGNITUDE DSPLV 

5450 PRINT "TYPE COMMAND PH, H6, RSJ, EX" 

5478 INPUT Ft ■* 

5480 IF FS * “pH" THEN GOTO 5558 

5490 IF FS * "M6“ THEN GOTO 5650 

5500 IF FS * "EX" THEN GOTO 5520 

5505 IF FS * -Rfci- THEN GOTO 5511 

5510 GOTO 5450 

5511 PRINT “WHICH SPECTRAL COMPONENT C 9-127 )?" 

5513 INPUT HS 

5514 IF HS < B OR HS > 127 THEN PRINT "ERROR: COMPONENT RANGE IS 8 TO 12 
7, INPUT AGAIN" • GOTO 5513 

5515 IF K S a “fiS" THEN PRINT* “REAL PART: "VRCMS): PRINT "IM. PART: "VICW 
S> 

5517 PRINT "REAL PART: "ZRCHS>: PRINT "IM. PART: "ZICWS) 

5519 GOTO 5458 

5529 RETURN 

5546 REH PHASE 

5550 HGR : HCOLOR* 7 

5560 HPLOT 0,0 TO 0,159 

5562 FOR V * 144 TO 0 STEP - 18 

5564 HPLOT 8,V TO 3,Y 

5566 NEXT V 

5568 HPLOT 0,159 TO 272,159 
5570 FOR X a is JO 272 STEP 8 
5572 HPLOT X,l53 TO X,156 

5574 NEXT X 

5575 GOSUB 5815 

5576 HPLOT 16,72 TO 272,72 

5578 FOR X = 8 TO 255 STEP 2 

5579 V = X / 2 

5588 Z a 72 - P<v> * 22.9 
5582 HPLOT X + 16,72 TO X + 16,2 
5584 NEXT X 

5612 IF KS * "ft$" THEN PRINT "****** PHASE PLOT: "RCHS" * ♦ ** »• “» GOTO 56 
18 

5614 IF E = 99 THEN PRINT "*** * «»» « POWER SPEC: "RCHS" SOTO 5 

618 

56 15 PRINT "** X-POWEP SPEC: "RCHS", "XCHS" **" 

5618 GOSUB 780 

5628 GOTO 5450 

5648 REH MAGNITUDE 

5650 HGR : HCQLQR* 7 

5669 HPLOT 8,0 TO 0,159 

5662 FOR V a 144 TO 8 STEP - 8 

5664 HPLOT 8,V TO 3,V 
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5666 NEXT V 

5668 HPLOT 0,159 TO 272,159 
5670 FOR X « 16 TO 272 STEP 8 
5672 HPLOT X,159 TO X,156 
5674 NEXT X 

5676 HPLOT 16,144 TO 272,144 
5678 60SUe 5815 

5630 PRINT "SUPRESS SPECTRAL COMPONENT SC V/N >?" 

5690 INPUT ES 

5700 IF E S = "V" THEN SOTO 5722 

5710 IF ES * -N- THEN SOTO 574f 

5720 60T0 5688 

5722 PRINT “WHICH COMPONENTC 0-127)?“ 

57’^4 I NP UT CP 

5725 IF CP < 0 OR CP > 127 THEN PRINT “ERROR: COMPONENT RAN6E IS 0 TO 12 

7, INPUT AGAIN*': SOTO 5724 

5726 W-. CP) * 0 
5728 SOTO 5680 

5741 PRINT “FULL, T4ALF, OR QUARTER DSPLVCF/H/Q >?“ 

5742 INPUT SOS 

5743 IF SOS * “F" THEN SO » 2s SOTO 5748 

5744 IF SOS * “H" THEN SO = 4: SOTO 5748 

5745 IF SOS * “Q“ THEN SO * 8: GOTO 5743 

5746 SOTO 5741 

5743 PRINT “SCALING FACTOR?" 

5749 INPUT A 

5750 FOR X * 0 TO 255 STEP SO 
5755 V = X / SO 

5760 Z * 144 - M<V) / A 
5770 IF Z > 144 OR 2 < 0 THEN GOTO 5748 
5790 HPLOT- X + 16,144 TO X + 16, Z 
5880 NEXT 

5805 IF ES * “V“ THEN H<0) * K 

5807 IF K S = “AS" THEN PRINT “****** HAS. PLOT: “RCHS“ ***■***": SOTO 561 
8 

5810 SOTO 5614 

5815 FOR V * 158 TO 156 ST EP - 1 

5820 HPLOT 1,V TO 16,V: NEXT 

5825 FOR X * 1 TO 3 

5830 HPLOT X,158 TO X,143 

5335 NEXT 

5840 RETURN 

5850 REM CORRELATION 

5860 PRINT “SUPRESS ANV SPECTRAL COMPONENTS^ V/N >?“ 

5865 INPUT SCS 

5870 IF SCS = “V THEN SOTO 5880 
5875 IF SCS * M N" THEN 60T0 5918 
5878 GOTO 5360 

5880 PRINT "WHICH COMPONENTC 1-128 >?“ 

5885 INPUT WC 

5837 IF WC < 1 OR WC > 128 THEN PRINT “ERROR: COMPONENT RANGE IS 1 TO 12 

8, INPUT AGAIN" : GOTO 5885 
5898 ZR< WC > = 0: ZIC WC > = 0 
5895 SOTO 5860 

5910 AR< 0 ) * 0:AI(0 > a 0 
5920 FOR X a i TO 128 
5930 V a 256 - X 
5940 ARC X > a 2P.CX) 

5950 ARC V) * ZRCX) 

5960 AIC V ) a - 2ICX) 

5970 AIC X ) a ZICX) 

5980 NEXT 
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5990 FRO » -1 
9000 80SUB 3000 
8605 H * 80 
6006 60SUB 1200 

6008 IF E « 99 THEN PRINT "****** ftUTO-CORR: 
6010 PRINT "*■#■* X-CORRs “RCHS“, "XCHS" *#*“ 
6812 PRINT “DSPLY CORR-FN flGflIN<YxN>?“ 

6014 INPUT PS 

6016 IF PS = "Y“ THEN SOTO 6806 
6018 IF PS « ~N“ THEN SOTO 6020 


8019 SOTO 6012 
6020 60SUB 700 
6036 SOTO 5850 


"RCHS“ #*+** * w ; SOTO 6012 
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F. Demonstration Run 


ORIGINAL PAGE IS 
OF POOR QUALITY 


This Appendix consists of a listing of an actual run of DAISE printed 
by the SILENTYPE* printer. The user is urged to read and follow the 
execution carefully, since it demonstrates vertually all the capabilities of 
DAISE and the routes leading there as well as the effects of certain deci- 
sions taken during the execution. The raw data sampled is a 10 Hz burst 
sinesoidal signal with half duty cycle. It was picked up by all microphone 
systems using Globe microphones, hence the file name. 

The sampling frequency constant was inputed as 15, implying 128 Hz fre- 
quency, that leaves us with the following resolutions: 

128 

Frequency resolution ■ — =— ■ .5 Hz 

256 

Time resolution ■ — ■ .0078 seconds. 

128 

Notes 

1 DAISE can not print any raw data and what is presented on the next page 
was printed by a different program for demonstration purposes. 

2 Whenever a scaling factor is inputed throughout the run, a HGR plot is 
displayed on the video screen and the message appearing after the scaling 
factor is what gets printed on the video screen underneath the display. 
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****RE«L OPTS: GL0B2-C1, QR-1 **** 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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I EXEC MASTER #01 


1 BftUN Q A I SE . OBJ #01 

Tv&£ COMMAND SR# PR# AN# C A# EX# FOR: 
SAMPLE /PREMIER# ANAL V$I$#CflTAL00#EXI7 

TSA 

INPUT SAMPLING FREQUENCY CONSTANT 
?15 

TYPE FILE NAME: 

?GL0E'2 

TYPE COMMAND SA# PR# AN# CA# EX# FOR: 
$AMPLE#PP£UIEX#ANALVSJ$#CATALGG#EXIT 

?CA 


ORIGINAL PAGE & 
OF POOR QUALITY 


DISK UOLUME 254 

A 002 HELLO 
B 006 SL032-C1 
B 606 6L0B2-C2 
B 006 GLGB2-C3 
B 006 6L0B2-C4 
B 006 6L0B2-C5 
B 006- GL0B2-C6 
B 006 8L0B2-C7 
B 006 SL0B2-C8 

TYPE COMMAND SA# PR, AN, CA# EX# FOR: 

$AMPL£#PREUI£H#ANALYSI$#CATAL06#EXIT 

7 Pft 

DESIRED FILE AND CHANNEL? 

7 8L0B2-C1 

DESIRED QUART ER< 1-4 >? 

?i 

SCALING FACTOR? 

n 

**** REAL DATA: GL0B2-C1# QR-1 **** 
TYPE COMMAND SA, PR# AN, CA# EX# FOR: 
SAMPLE# PREU I EH# ANALYSIS# CATALOG# EX IT 
?PF: 

DESIRED FILE AND CHANNEL 7 
?6lOB2-C2 

DESIRED QUARTER< 1-4 >? 

?1 

SCALING FACTOR? 


REAL DATA: GL0B2-C2 # QR-1 
t VPE CCtMMPND SA, PR, AN# CA# EX# "OR: 
SAMPLE # PREU I EH # ANAL YS I S # CATALOG # EX I T 

7 AN 

TYPE COMMAND SE# CO# DI# EX, FOR: 
SET-UP, COMPUTE# DISPLAY# EXIT 

7 SE 

INPUT REFRENCE CHANNEL 
7 GL 062-Cl 

TYPE COMMAND AC# CC# FOP: 
AUTO-CORRELATION # CROSS-CORRELAT I ON 

?CC 

INPUT X -CHANNEL 
■'GL0B2-C2 

STARTING BLOCKS. 0-7 >? 

*o 

HOW many BLOCKSC 1-8 .‘ 7 


* 

ji 
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TYPE COMMAND SE, CO, DI, EK, FOR: 

SET-UP, COMPUTE# DISPLAY, EXIT 

rco 

TYPE COMMAND SE, CO, 01, EX, FOR: 

SET-UP, COMPUTE* DISPLAY, EXIT 
TO I 

TYPE- COMMAND PS, AS, CR, EX, FOR: 

POWER SPEC, AMPL. SPEC, CORRELATION, EXIT 

?-AS 

TYPE COMMAND PH, MS, R&I, EX 

Tpw 

phase plot: slgbs-ci » »**» ♦ 

MOULD YOU LIKE TO PRINT THE DISPLAY< Y/N >? 
?Y 


PHASE OF AMPL. SPEC***GL0B2-CHhh* AUG 



TYPE COMMAND PH, MS, Ril, EX 

TM© 

SUPRESS SPECTRAL COMPONENTS*: YrN >? 

?v 

WHICH COMPONENT < 0-127 >? 

?0 

SUPRESS SPECTRAL COMPONENTS Y^N )? 

?N 

FULL, HALF, OR QUARTER OSPLYCF^H/Q)? 

TP 

SCALING FACTOR? 

71 

MAS. PLOT: SLQB2-C1 » « ■*•» » * 

MOULD YOU LIKE TO PRINT THE DISPLAY< V.'NT? 
TV 
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AMPL. SPEC . ***6L0BZ.'-C 1 ***4 AUG***SCALE*l 


GkIU^AL PAG!: 19 

OF POOR QUALITY 



TYPE CQMMANO PH, MG, Rtl, EX 
?MS 

SUPRESS SPECTRAL COMPONENTS* V/N )? 

?N 

FULL, HALF, OR QUARTER DSPLV< F/H/'Q >? 

?H 

SCALING FACTOR? 

?. 5 

SCALING FACTOR? 

****** MAG. PLOT: 6L0B2-C1 *** »** 

MOULD YOU LIKE TO PRINT THE DISPLAY* V/hO? 
■N 

TYPE COMMAND PH, MG, RU , EX 
VMS 

SUPRESS SPECTRAL COMPONENTS* Y/N >? 

?N 

FULL, HALF, OR QUARTER 0SPLY< F/H/Q P 
?H 

SCALING FACTOR? 

•.A 
/ r 

»**«*■* MAG. PLOT: 6L0S2-C1 * x x *** 

MOULD YOU LIKE TO PRINT THE DISPLAY* V/N >* 

7V 
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1 « « ‘ » 1 I - t L i. -i-1— 1— I— 1— L-JL JL 


i 


AMPL. SPEC . *-**6L0B2-C 1 AUG***SCALE=2 


CRiGIKAL & 

OF POOR QUALITY 




L 


ii 


J ! . ( ■ 1 n ■ 1 . « 1 1 


tl , « 


■ » t «--» * lit i . 


J— L 


_i l i i 


TYPE COHHANO PH, H6, R&I, EX 

TEX 

TYPE COMMAND PS, AS, CP, EX, FOR: 

POWER SPEC, AMPL. SPEC, CORRELATION, EKIT 
TPS 

TYPE COMMAND PH, MS, R&I, E* 

tmr 

SUPRESS SPECTRAL COMPONENTS* Y^N >? 

?y 

MH I CH COMPONENTS 0- 1 2? >? 

T0 

SUPRESS SPECTRAL COMPONENTS': Y-'N )? 

?N 

FULL, HALF, OR QUARTER OSPLY* F'HsQ >? 

7 H 

SCALING FACTOR? 

“'10 

SCALING FACTOR? 

"100 

** X -POWER SPEC: GL0B2-C1, GL0B2-C2 ** 
WOULD YOU LIKE TO PRINT THE DISPLAYS Y/N>? 
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X-PQMER***GL0B2-C 1 CROSS GL0B2-C2***4 AUG***SCALE=100 


r 

[ 

F 

r 

r 

r 

r 

h 

r 


« ■ « 






rajsE 

0F POOR QUALITY 


■ - 1 1 1 1 1 


TVPE COMMAND PH, MG, RSJ, EX 
?H6 

SUPRESS SPECTRAL COMPONENTS*! Y/N )? 

?Y 

WHICH COMPONENT* 8-12? >? 

?1 

SUPRESS SPECTRAL COMPONENTS* Y/N >? 

?N 

PULL, HALF, OR QUARTER DSPLWF/H/Q)? 

SCALING FACTOR? 

?10 

SCALING FACTOR? 

** X-POWER SPEC: GLGB2-C1, SLG62-C2 ** 
WOULD YOU LIKE TO PRINT THE DISPLAY*. V/N )? 

?N 

TVPE COMMAND PH, MG, Rfcl, EX 
?MG 

SUPRESS SPECTRAL COMPONENTS* V/N >? 

?N 

FULL, HALF, OR QUARTER DSPLY* F/H/Q )? 

?Q 

SCALING FACTOR? 

: '20 

** X— POWER SPEC: 6LGB2-C1 , GL0B2-C2 
WOULD YOU LIKE TO PRINT THE D I SPLAY* Y/N )? 
"V 
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X-P0HER***GLQB2-Cl CROSS GLOB2-C2***4 AU6***SCAL£=20 


r 

r 


OWGfMn r* la 
OF POOR QUALITY 



I ■ 1 t - t L 

t I I 1 1. I U. 1-1 -k 


k + l-i-1 

J l I I L 




i 


i I I L 


0mM 

J L 


TYPE COMMAND PH# HG# R&I# EX 

mi 

HHICH SPECTRAL COMPONENT-: 0-12? )? 

?1 

REAL PART: 6421.46912 
IM. PART: 2501.65725 ' 

TYPE COHMANO PH# M6# R*I# EX 

mi 

WHICH SPECTRAL COMPONENT-: 0*127 )? 

720 

REAL PART: -1292.40875 
IM. PORT: -2991.91292 
TYPE COMMAND PH, M6# R&I# EX 

7EX 

TYPE COMMAND PS# AS# CR, &i, FOR: 

POWER SPEC# AMPL. SPEC# CORRELATION# E^IT 

TCP 

SUPRESS ANY SPECTRAL COMPONENTS*: Y/N >? 

?N 

SCALING FACTOR? 

•T 

SCALING FACTOR? 

•?4 

SCALING FACTOR? 

?19 

*** X-CORR: GL0B2-C1# GL0B2-C2 **# 

OSPLV CORR-FN AGAIN* V/N>? 

7 V 

SCALING FACTOR? 

? 1 9 

•hh* x-CORR: GL0B2-C1 # GL0B2-C2 *-** 

OSPLV CORR-FN AGAIN*Y/N V? 

7N 

WOULD YOU LIKE TO PRINT THE DISPLAY* Y/N>? 



OKiiafft.'X P/',-"'-:’ tg 

OF POOR QUALITY 


X-C0RR***6L0B2-C1 CROSS GL0B2-C2***4 AUG*-**SCALE*10 



TVPE COMMAND PS, AS, Cfi, EX, FOR: 

POHER SPEC, RMPL. SPEC, CORRELATION, EXIT 
?CR 

SUPRESS ANY SPECTRAL COMPONENTS 1 ' Y.'N )? 

AH I iX COMPONENT •' 1-123 >? 

SUPRESS ANY SPECTRAL COMPONENTS^. Y/N >? 

?N 

SCALING FACTOR? 

•-5 

~ X-CORR: 6L0B2-C1, GL0B2-C2 *** 

DSPLY CORR-FN A6AIN< Y/N)? 

■•N 

MOULD YOU LIKE TO PRINT THE DISPLAY< V/N >? 
?Y 


X-C0RR***6LQB2-C1 CROSS SL0B2-C2***4 AUS***SCALE=5 



TYPE COMMAND PS# AS, CR# EX# FOR: 

POWER SPEC# AHPL. SPEC# CORRELATION# EXIT 

?CR 

SUPRESS ANY SPECTRAL COMPONENTSC V^N >? 

?V 

WHICH COMPONENT < 1 - 1 28 >? 

?17 

SUPRESS ANY SPECTRAL COMPONENTS< Y/N '•? 

?Y 

WHICH COMPONENTC 1-128 >? 

?18 

SUPRESS ANY SPECTRAL COMPONENTSC V/N )? 

7 Y 

WH I CH COMPONENT C 1-123 >? 

?19 

SUPRESS ANY SPECTRAL COMPONENTSC Y/N >? 

~Y 

WH I CH COMPONENT < 1-123 >? 

?21 

SUPRESS ANY SPECTRAL COMPONENTSC Y/N >? 

?Y 

WHICH COMPONENT C 1-1 23 >? 

?22 

SUPRESS ANY SPECTRAL COMPONENTSC Y/N Y< 

?Y 

WHICH COMPONENT C 1-123)? 

SUPRESS ANY SPECTRAL COMPONENTSC Y/N >? 

?N 


ORIGIN! PASS K 

SCALING FACTOR? OF POOR QUALITY 

?5 

*** X-CORR: GLGB2-C1, 6L0B2-C2 *** 

OSPLV CORR-FN AGAINCV/N)? 

?N 

WOULD VOU LIKE TO PRINT THE DISPLAY Y/N>? 

?Y 


X-C0RR‘5hhc6L0B2-C1 CROSS GL0B2-C2***4 AUG***SCALE=5 



TYPE COMMAND PS* AS, CR, EK, FOR: 

POWER SPEC, AMPL. SPEC, CORRELATION, EXIT 
TPS 

TYPE COMMAND PH, M6, R&I, EX 
TM6 

SUPRESS SPECTRAL COMPONENTS*: Y/N >? 

?N 

FULL, HALF, OR QUARTER DSPLYC F/H/Q >? 

?H 


oki-cHr;.. 

OF POOR QUALITY 


SCALING FACTOR? 

?20 

SCALING FACTOR? 

?*3 

SCALING FACTOR? 

?180 

SCALING FACTOR? 

?1000 

SCALING FACTOR? 

?4000 

SCALING FACTOR? 

?10000 "« 

** X-PGWER SPEC: 6LQB2-C1, GL0B2-C2 ** 
WOULD YOU LIKE TO PRINT THE DISPLAY* Y/N y? 

TV 


X-P0WER#*#6L0B2-C1 CROSS SLOB2-C2***4 AUG***SCALE=10000 


F 


F 

r 
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TYPE COMMAND PH, MG. RU. EX ORlUfK/iL P , : ~. s , ^ 

SUPRESS SPECTRAL COMPONENTS* Y/N >? 0F P °°* QUALITY 

?Y 

WHICH COMPONENT* 0-12? )? 

SUPRESS SPECTRAL COMPONENTS* Y/N )? 

FULL, HALF, OR QUARTER OSPLY* F/H/Q >? 

?H 

SCALIN6 FACTOR? 


?20 

•5Hr X-POWER SPEC: 6L0B2-C1, GL0B2-C2 ** 
WOULD YOU LIKE TO PRINT THE DISPLAY* Y/N >? 


?Y 


X-POWER*hhj6LOB2-C 1 CROSS 6LDB2-C2* 1 *4 AUC* * kSCALE-20 



TYPE COMMAND PH, 

TEX 

TYPE COMMAND PS, 
POWER SPEC, AMPL. 


MG , R& l , E» f: i 

AS, CR, EX, FOR: 

SPEC, CORRELATION, EXIT 


TEX 

TYPE COMMAND SE, 
SET-UP, COMPUTE, 


CO, 01, EX, FOR: 
DISPLAY, EXIT 


! t A 


TYPE COMMAND SA, PR, AN, CA, EX, FOR: 
SAMPLE , PREU I EH , ANALYS IS, CATALOG .EXIT 



